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PARTIALLY PENALIZED IMMERSED FINITE ELEMENT 
METHODS FOR ELLIPTIC INTERFACE PROBLEMS * 

TAO LINt, YANPING LIN*, AND XU ZHANG§ 


Abstract. This article presents new immersed finite element (IFE) methods for solving the 
popular second order elliptic interface problems on structured Cartesian meshes even if the involved 
interfaces have nontrivial geometries. These IFE methods contain extra stabilization terms intro¬ 
duced only at interface edges for penalizing the discontinuity in IFE functions. With the enhanced 
stability due to the added penalty, not only these IFE methods can be proven to have the optimal 
convergence rate in an energy norm provided that the exact solution has sufficient regularity, but 
also numerical results indicate that their convergence rates in both the fJ^-norm and the L^-norm 
do not deteriorate when the mesh becomes finer which is a shortcoming of the classic IFE methods 
in some situations. Trace inequalities are established for both linear and bilinear IFE functions that 
are not only critical for the error analysis of these new IFE methods, but also are of a great potential 
to be useful in error analysis for other related IFE methods. 
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1. Introduction. Without loss of generality, we consider a domain 17 that is a 
union of rectangular domains in and we assume that 17 is formed by two different 
materials separated by a curve F. In particular, this means F separates 17 into two 
sub-domains 17“ and 17+ such that 17 = 17“ U 17+ U F. Consequently, the diffusion 
coefficient /3 on 17 is assumed to be a piecewise constant function: 


P{x,y) 


13 , (a;, 2 /) G 17 , 
13+, (a:,y)Gl7+, 


such that min{/3“,/3+} > 0. 

The main purpose of this article is to present a group of partially penalized im¬ 
mersed finite element (IFE) methods using Cartesian meshes to solve popular elliptic 
interface problems appearing in many applications in the following form: 

(1.1) -V ■ {l3Vu{x,y)) = f{x,y), (a:, y) G 17“ U 17+, 

(1.2) u{x,y)=0, {x,y)€dn, 


together with the jump conditions on the interface F C 17: 


(1.3) 

(1.4) 


M |r = 0, 

,9it1 

1 "=“' 


The homogeneous boundary condition (1.2) is discussed here for simplicity’s sake, the 
method and related analysis can be readily extended to interface problems with a 
non-homogeneous boundary condition. 
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A large number of numerical methods based on Cartesian meshes have been intro¬ 
duced for elliptic interface problems. Since Peskin’s pioneering work of the immersed 
boundary method [49], a variety of methods have been developed in hnite difference 
formulation, such as the immersed interface method [35], the matched interface and 
boundary method [61], and the ghost fluid method [24]. We refer to the book [37] for 
an overview of different numerical methods in finite difference framework. 

In finite element formulation, certain types of modifications need to be executed 
for elements around the interface. One way is to modify the weak formulation of finite 
element equations near the interface. We refer to some representative methods such 
as the penalty hnite element method [4, 11], the unhtted hnite element method [28], 
the discontinuous Galerkin formulation methods [12, 27]. An alternative approach is 
to modify the approximating functions around the interface, for instance, the general 
hnite element method [6, 7], the multi-scale hnite element method [18, 23], the ex¬ 
tended hnite element method [21], the partition of unity method [8, 9], to name just 
a few. 

Immersed hnite element (IFE) methods are a particular class of hnite element 
(FE) methods belonging to the second approach mentioned above, and they can solve 
interface problems with meshes independent of the interface [1, 2, 29, 30, 32, 33, 36, 
39, 42, 47, 54, 55]. If desired, an lEE method can use a Cartesian mesh to solve a 
boundary value problem (BVP) whose coefficient is discontinuous across a curve F 
with a non-trivial geometry. The basic idea of an IFE method is to employ standard 
FE functions in non-interface elements not intersecting with the interface F, but on 
each interface element, it uses IFE functions constructed with piecewise polynomials 
based on the natural partition of this element formed by the interface and the jump 
conditions required by the interface problem. The lEE functions are macro-elements 
[13, 20], and each lEE function partially solves the related interface problem because 
it satishes the interface jump conditions in a certain sense. Also, the lEE space 
on an interface element is consistent with the corresponding EE space based on the 
same polynomial space in the sense that the lEE space becomes the FE space if the 
discontinuity in the coefficient j3 disappears in that element, see [29, 30] for more 
details. 

IFE methods have been developed for solving interface problems involving sev¬ 
eral important types of partial differential equations, such as the second order elliptic 
equation [1, 2, 25, 29, 30, 33, 34, 36, 39, 41, 54, 55, 59, 60], the bi-harmonic and beam 
equations [42], the planar elasticity system [26, 40, 46, 47], the parabolic equation 
with hxed interfaces [3, 45, 56], and the parabolic equation with a moving interface 
[32, 43, 44]. When jump conditions are suitably employed in the construction of lEE 
functions for an interface problem, the resulting lEE space usually has the optimal 
approximation capability from the point view of polynomials used in this IFE space 
[1, 15, 16, 30, 38, 52, 60]. Numerical examples [2, 36, 38, 39] demonstrate that methods 
based on lEE spaces can converge optimally for second order elliptic interface prob¬ 
lems. However, the proof for their optimal error bounds is still elusive except for the 
one dimensional case [2], even though there have been a few attempts [17, 31, 34, 57]. 
Eor two dimensional elliptic interface problems, only a suboptimal convergence in the 
iJ^-norm has been rigorously proven [31]. 

One of the major obstacles is the error estimation on edges between two interface 
elements where lEE functions have discontinuity. Certain types of trace inequalities 
are needed and can be established, but it is not clear whether the generic constant 
factor in these inequalities is actually independent of the interface location. The 
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scaling argument in the standard finite element error estimation is not applicable 
here because the local IFE spaces on two different interface elements are not affine 
equivalent in general. Besides, numerical experiments have demonstrated that the 
classic IFE methods in the literature often have a much lager point-wise error over 
interface elements which, we believe, is caused by the inter-element discontinuity 
of IFE functions. In some cases, the convergence rates can even deteriorate when 
the mesh becomes finer. These observations motivate us to apply a certain penalty 
over interface edges for controlling negative impacts from this discontinuity. Natural 
candidates are those well known penalty strategies for handling inter-element dis¬ 
continuity in interior penalty Galerkin methods and discontinuous Galerkin methods 
[5, 10, 14, 22, 48, 51, 53, 58]. These considerations lead to the partially penalized IFE 
methods in this article. Theoretically, thanks to the enhanced stability by the penalty 
terms, we are able to prove that these new IFE methods do converge optimally in an 
energy norm. In addition, we have observed through abundant numerical experiments 
that these partially penalized IFE methods maintain their expected convergence rate 
in both Fl^-norm and L^-norm when their mesh becomes finer and finer while the 
classic IFE methods cannot maintain in some situations. 

The partial penalty idea has also been used in the unfitted finite element method 
[28]. In this method, penalty terms are introduced on interface instead of interface 
edges because approximating functions are allowed to be discontinuous inside interface 
elements but they are continuous on element boundaries within each subdomain. 
IFE methods reverse this idea by imposing continuity of approximating functions 
inside each element but allowing discontinuity possibly only across interface edges. In 
addition, on the same mesh, the unfitted finite element method has a slightly larger 
number of degrees of freedom than IFE methods. On the other hand, the unfitted 
finite element method has been proven to have the optimal convergence rate under the 
usual piecewise regularity [28] while the analysis in the represent article needs to 
assume a piecewise or regular in order to establish the optimal convergence 

for the partially penalized IFE methods. 

Also, we note that these partially penalized IFE methods and their related error 
analysis can be readily modified to obtain IFE methods based on the discontinuous 
Galerkin formulation with advantages such as adaptivity even with Cartesian meshes. 
However, on the same mesh, the DG IFE methods generally have far more global de¬ 
grees of freedom. For instance, on a Cartesian triangular mesh, a DG IFE method has 
about 6 times more unknowns than the classic IFE method. The partially penalized 
IFE methods presented here have the same global degrees of freedom as their classic 
counterparts; hence they can be more competitive in applications where advantages 
of DG IFE methods are not needed. 

The rest of this article is organized as follows. In Section 2, we derive partially 
penalized IFE methods based on either linear or bilinear IFE functions for the interface 
problem. In Section 3, we show that the well-known trace inequalities on an element 
are also valid for linear and bilinear IFE functions even though they are not 
functions locally in an interface element. In Section 4, we show that these IFE schemes 
do have the optimal convergence rate in an energy norm. In Section 5, we will present 
numerical examples to demonstrate features of these IFE methods. 

2. Partially penalized IFE methods. Let Th, 0 < /i < 1, be a family of 
Cartesian triangular or rectangular meshes on H. For each mesh Th: we let J\fh be the 
set of vertices of its elements, and let Eh be the set of its edges and Eh be the set of 
interior edges. In addition, we let Tf^ be the set of interface elements of Th and 7^" be 
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the set of non-interface elements. Similarly, we let be the set of interior interface 
edges and let be the set of interior non-interface edges. For every interior edge 
S G f/i, we denote two elements that share the common edge B by Tb,i and Tb, 2 - 
For a function u defined on Tb.i U Ts 2 , we denote its average and jump on B by 

{u}b = -^{{u\tb,i)\b + {u\tb,2)\b), [u]b = (mItb.JIb “ {u\tb,2)\b- 

For simplicity’s sake, we will often drop the subscript B from these notations if there 
is no danger to cause any confusions. We will also use the following function spaces: 

= {v G I u|n» G W’’^(0®), s = -I- or —} for r > 1 and 1 < p < 00, 

equipped the norm 

+ lkll^^.p(n+), Vu G 

As usual, for p = 2, we use H^{n) = and denote its corresponding norm by 

ll^llr = IMlI-cn) = ^ 

With a suitable assumption about the regularity of F and / (e.g. [4]), we can 
assume that the exact solution u to the interface problem is in To derive a 

weak form of interface problem described by (1.1)-(1.4) for an IFF method, we will 
use the following space: 

Vh = {v \ V satisfies conditions (HV1)-(HV4) described as follows} 

(HVl) v\kGH\K), 'iKGTh. 

(HV2) V is continuous at every X G Aft- 
(HV3) V is continuous across each B G 
(HV4) r;|9o = 0. 

We multiply equation (1.1) by a test function v G Vh, integrate both sides on each 
element K G Th, and apply Green’s formula to have 


/3Vu • VudX — / /3Vm • nvds = / vfdX. 


IK 


IdK 


'K 


Summarizing over all elements leads to 

(2.1) 'S^ I jdXv-XudX— j • ns} [u]ds = f vfdX. 


Beei 


Here we have used the fact that 

{/3Vm • nslg = (/3Vu • nB)|B, VH G 

Because of the regularity of u, for arbitrary parameters e, a > 0, and > 0, we have 


( 2 . 2 ) 


e X] / {/3Vu • ns} [njds = 0, X 


Best 


Beei 


\B\ 


[n][u]ds = 0. 
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Therefore, adding (2.2) to (2.1) leads to the following weak form of the interface 
problem (1.1)-(1.4): 


(2.3) V / pVvVudX- V / {/SVu-nsIHds 


Ken 


Best 


+e V / {pVvnB}[u]ds+ V [ -^a[v][u]ds = [ vfdX, WveVh. 

Jii JB „ Jb \^\ JQ 

BeSh BeCy^ 

We now recall the linear and bilinear IFE spaces to be used in our partially penalized 
IFE methods based on the weak form (2.3). On each element K G Th, we let 


Sh{K) = span{cj}j{X), 1< j < die}, dK 


3, if A" is a triangular element, 

4, if AT is a rectangular element. 


where , 1 < j < dx are the standard linear or bilinear nodal basis functions for 
K G 7)j"; otherwise, for K G 7^, < j < dx are the linear or bilinear IFE basis 

functions discussed in [38, 39] and [30, 41], respectively. Then, we define the lEE 
space over the whole solution domain as follows: 


Ship) = {n 1 n satisfies conditions (lEEl) - (IFE3) given below} 

(IFEl) v\xGShiK), WKGTh. 

(IFE2) V is continuous at every X G Afh- 
(IFE3) tijan = 0. 

It is easy to see that Shp) C 14(0). Now, we describe the partially penalized 
IFE methods for the interface problem (1.1)-(1.4): find Uh G Shp) such that 

(2.4) ahivh,Uh) = (vh, f), VvhGShp), 

where the bilinear form a^(-, •) is defined on Shp) by 


ahivh,Wh)= 'Y' / jdVvh-VwhdX - 'Y / {jdVWh ■ ns} [vh]ds 


Ken 


Beet 


(2.5) + e ^ / {PVvh ■ ns) [wh]ds + Y 


Besi 


Best 


\B\ 


dvh][wh]ds, 'dvhyWh G Shp). 


3. Trace inequalities for IFE functions. Using the standard scaling argu¬ 
ment, we can obtain the following well known trace inequalities [50]: there exists a 
constant C such that 

(3.1) < C \bY (lkllL2(^) + h 11 , Va £ H\K), 

(3.2) IIVz;ll^.(5) < C \kY^ (II+ h ||, Va G H^{K). 
where B is an edge of AT. 

Our goal in this section is to extend these trace inequalities to IFE functions in 
Sh{K) for AT G Th- Eirst, we recall that ShiK) C C{K) 0 iJ^(Ar) for all K G Th 
[30, 38]. This implies that inequality (3.1) is also valid for v G ShiK) even if AT G 7)}. 
However, the second trace inequality (3.2) cannot be applied to G ShiK) with 
K GTh because v ^ iJ^(Ar) in general. 
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3.1. Trace inequalities for linear IFE functions. It is relatively easier to 
prove that the trace inequality for a linear IFE function in a triangular interface 
element is true because its gradient is a piecewise constant function. Without loss of 
generality, we consider the following triangular interface element 

K = AA 1 A 2 A 3 , = (0,0), A 2 = (h, 0), A 3 = (0, h). 


Assume that the interface F intersects the edge of K at points D and E and the 
straight line DE separates K into K~ and Ar+, see the illustration on the left in Fig. 
3.1. Consider a linear IFE function on K in the following form 


(3.3) 


v{x,y) 


V {x,y) = c-^ + c^x + c^y, ii {x,y) £ K , 
v+{x,y) =c'^ + c^x + cjy, A{x,y)eK+, 


which satisfies the following jump conditions [38]: 


(3.4) 


-{D)=v+{D), v-{E)=v+{E), p- 


dv 


= /3+ 


dv+ 




Fig. 3.1. a triangular interface element (left) and a rectangular interface element (right) 

Lemma 3.1. There exists a eonstant C > 1 independent of the interface location 
such that for every linear IFE function v on the interface element K defined in (3.3) 
the following inequalities hold: 


(3.5) 


C 


ici,4,cf)\\ < ||(Ci ,C 2 , 03 )]] < C\\{cf,cJ,cJ) 


Proof. We prove the second inequality in (3.5), and similar arguments can be 
used to show the first one. Applying the jump conditions (3.4) we can show that 
coefficients of v{x,y) must satisfy the following equality: 



where A/®, s = —, + are two matrices. Without loss of generality, we further assume 


D = (0, d/i), E = {eh, 0), with 0<(i<l,0<e<l. 
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Then 


/l 0 dh \ 

M® = jl eh 0 j,s = — or + 

\0 —j3^dh —13‘^ehJ 

whose determinant is 

det(M®) = + e'^)h^, s = - or + 

which is nonzero because (d, e) 7 ^ (0, 0). Hence we can solve for cf s in terms of c~ to 
have 


(3.6) 

(3.7) 
with 


— fil^l “t“ fi‘1^2 ^ — 1; 2, 3, 

_ 9ijP ^ . g 

/3+(d2 + e2) + ^+(d2+e2)’ 


dll = 0 , 5 n = + e^, = -d^eh.g^^ = d^eh, = -de^h,g^^ = de^h, 

921 ~ dji ~ 0, (^22 ~ d , g 22 — e , 523 “ de, g^^ = —de, 

9zi — Oj dai ~ 932 ~ de,g ^2 — ~de, 533 = e ,( 7 ;^ = d . 

Therefore, there exists a constant C that depends on /3“ and /3+, but is independent 
of d, e, such that 


\h\ <C, 1< i,j < 3. 

Then, the second inequality in (3.5) follows from (3.6) and the above bounds for 
1/^1, l<i,J<3. □ 

Remark 3.1. In the proof of Lemma 3.1, we have shown fn = 0,i = 2,3. 
Consequently, we can show that there exists a constant C > 1 such that 

(3.8) ^11(4:4)11 < ||(c 2 -,C 3 -)|| <C||(c+,c+)||. 

Now, we establish the trace inequality on a triangular interface element K = AH 1 H 2 H. 3 . 


Lemma 3.2. There exists a constant C independent of the interface location such 
that for every linear IFE function v on K the following inequalities hold 


(3.9) 

(3.10) 




L^K) 




, P = x,y, 


LHK) 


Proof. Without loss of generality, we assume again that the interface T intersects 
with the boundary of K at 

D = (0, dh), E = {eh, 0), with 0<d<l,0<e<l, 

and the line DE separates K into two subelement K~ and K'^ with A 3 G K~^ 
and ^ \K\. Furthermore, we assume B = A 1 A 3 is an interface edge with 
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B = B~ U B^. Similar arguments can be applied to establish the trace inequality in 
other cases. 

By direct calculations, we have 


11/3^ 


Il 2 (b+) 


= {ctf \B+\ (1/3+1 < {{ctr + (c+f) \K+\ 1/3+1) 


— R+ 


\B+\ 

\K+\ 


\/^Vz 


L^{K+) \K+\ 


V^V 


\K+\ 

2 


LHK) 


ie., 

(3.11) 






LHK) ■ 


Similarly, we can show that 


(3.12) 




x/^Vu 


L^{K) ■ 


On B , applying the estimates in Remark 3.1, we have 


/3 Vx 


2 

L^B-) 


{c^r \B-\ ( \r\f < C{{4f + (ctf) \K+\ |/3-| 


0/3- 


| 11 -| 

\K+\ 


x/ /3+Vu 


2 

L 2 (ix+) 


< 


c 


lll-l 

li^+l 


V^Vu 


2 

L 2 (if) ■ 


Hence, 


(3.13) 
Similarly, 

(3.14) 


/3-^y|L.(s_) <2Chi/2|x|-i/2 


x/^Vu 


L=(iX) ■ 


x/^Vu 


L 2 (/X) ■ 


Then, the combination of (3.11) and (3.13) yields the inequality (3.9) for p = x. 
Similarly, the inequality (3.9) for p = y follows from combining (3.12) and (3.14). 
Finally, (3.10) follows directly from (3.9). □ 


3.2. Trace inequalities for bilinear IFE functions. Without loss of gener¬ 
ality, we consider a rectangular interface element K with the following vertices: 


Hi = (0,0), H2 = (h,0), H3 = (0,h), A4 = {h,h). 

Again, assume that the interface F intersects with dK at points D and E and the 
linear DE separates K into two subelements K~ and A'+, see the illustration on the 
right in Fig. 3.1. We assume that K is one of the two types of rectangular interface 
elements [29, 30]: 

Type I interface element: The interface F intersects dK at 


D = (0, dh), E = {eh, 0), 0 < d < 1, 0 < e < 1. 
Type II interface element: The interface F intersects dK at 
D = {dh, h), E = {eh, 0), 0 < d < 1, 0 < e < 1. 
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On this interface element K, let n be a bilinear IFE function in the following form: 


(3.15) 


v(x,y) 


V {x,y) = + c^x + c^y + c^xy, ii {x,y) £ K , 
v+{x,y) = cf +c'^x + c^y + C 4 xy, if {x,y) e K+, 


which satisfies pertinent interface jump conditions [29, 31]. First, using similar argu¬ 
ments, we can show that the coefficients of a bilinear IFE function v satisfy inequalities 
similar to those in Lemma 3.1. 

Lemma 3.3. There exists a constant C > 1 independent of the interface location 
such that for every bilinear IFE function v on the interface element K defined in 

(3.15) the following inequalities hold: 

(3.16) ^ ||(c]^,cJ,c;]',C 4 )|| < ||(cr,C2 ,C;]',C4)|| < C||(c]^,cJ,C;]',C4)|| . 

The proof of trace inequalities for a bilinear IFE function is a little more complicated 
because its gradient is not a constant. The following lemma provides an aid. 

Lemma 3.4. Assume K is an interface element such that 

\in>l\K\, 


with s = — or -|-. Then there exists a polygon K C /F® and two positive constants Ci 
and C 2 independent of the interface location such that 


(3.17) 

(3.18) 


K 


> Cl \K\ 


L^k) 


>C2P^ihiC2f + h{clf + h^{c4f). 


Proof. Let us partition K into 4 congruent squares = 1,2,3,4 by the lines 
connecting the two pairs of opposite mid points of edges of K such that Ai is a vertex 
of Ki. Since 5 |7F|, one of these 4 small squares must be inside K^. Without 

loss generality, we assume that K 4 C . By direct calculations we have 


Ilw|li2(^^) = ^(12(c^)^-bl8c^c4/i-b7c^h2) > ^ 
ll^ylli2(A4) = ^(12 (c3)^ + 18c^C4/i -b Iclh"^) > ^ 


(^12 (c*)2 + (7-9ai)c2/i2 

^12 - (Cg)^ + (7 - QcT 2 )clh^ 


where ui and 172 are arbitrary positive constants. Letting ai = a £ (9/12, 7/9), i = 1, 2 
in the above inequalities leads to 

\W\l.^K.^>Ch\{c?2f + {clf+clh^) 

where 

C = min{12 - -, 2(7 - 9cr)} > 0. 
a 

Then, (3.17) and (3.18) follow by letting K = K 4 . □ 

Now, we are ready to establish the trace inequality for bilinear IFE functions on 
an interface element K = □^ 1 ^ 42243 ^ 4 . 
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Lemma 3.5. There exists a constant C independent of the interface location such 
that for every bilinear IFE function v{x,y) on K the following inequalities hold 


(3.19) 

(3.20) 


WPvph^B) < 


L^{K) 




, P = x,y, 


LHK) 


Proof. Without loss of generality, we assume that K is a. Type I interface element 
and B = A 1 A 3 is an interface edge. To be more specific, we also assume that A 4 G A"+ 
and |Ar“'"| > ^ |Ar|. Then 


B = B-U B+, B- = AiD, B+ = DA 3 . 

Direct calculations lead to 

|2 1 


(3.21) ||/3 = (/3 )‘^{dh{c2y + d'^h'^C2C4 +-d^h^cf) 

(3.22) ||/3 11^2(5-) = (i® ) dh{c3) , 

(3.23) ||/3+^;.||'2(b+) = (^+) 

(3.24) ||/3+^;y 11^2(5+) = (/5+)'(l - d)h{ctr. 

Applying (3.17) and (3.18) to (3.23) and (3.24) yields 

h 


(1 — d)h{c2)‘^ + (1 ~ (i^)/i^(cj)c4 + -(1 — d^)h^c\ 

o 


(3.25) 


W, II" <C- 


\P T 


P\\L^{B + ) - ^ 1^1 


h ^ 2 


< C- 

L^K+) - \K\ 


V^Vr 


LHK) 


, P = x,y. 


Moreover, applying (3.16), (3.17), and (3.18) to (3.21) and (3.22) leads to 


(3.26) V^V 


< 

LHK+) \K\ 




LHK) 


, p = x,y. 


Then, (3.19) follows from combining (3.25) and (3.26) together. Finally, (3.20) obvi¬ 
ously follows from (3.19). □ 

4. Error Estimation for Partially Penalized IFE Methods. We show that 
the IFE solution to the interface problem solved from (2.4) has an optimal convergence 
from the point of the polynomials used in the involved IFE spaces. Unless otherwise 
specified, we always assume that 7)i,0<ft.<lisa family of regular Cartesian 
triangular or rectangular meshes [19]. 

We start from proving the coercivity of the bilinear form ah{-,-) defined in (2.5) 
on the IFE space Sh{^) with respect to the following energy norm: 


(4.1) ||n^|U= f ^ / (3WvH-^vndX+ ^ [ 


yKen 


B 


^ [vh][vh]ds 


Besi 


) 


Lemma 4.1. There exists a constant k > 0 such that 


(4.2) 


K\\vh\\h < ah{vh,Vh), ^Vh G Sh{^) 
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is true for e = 1 unconditionally and is true for e = 0 or e = —1 under the condition 
that the stabilization parameter in is large enough. 


Proof. First, for e = 1, we note that the coercivity follows directly from the 
definitions of •) and IHI^- 
For e = —1, 0, note that 


(4.3) ah{vh,Vh)= ^ [ pVvh-yvhdX + [ {fiyvh ■ ns} [vh]ds 


Beet 


Beit 


\B\ 


r[vh][vh]ds, 


and the main concern is the second term on the right hand side. For each interface 
edge B G we let Kb^i € Th,i = 1, 2 be the two elements sharing B as their common 
edge. Then, by the trace inequality (3.10) or (3.20) and using a > 1, we have. 


/ {l 3 Vvh ■ ns} [vh]ds < UPVvh ■ nB}||L2(B) |IK]|Il2(b) 

J B 

< ^2■ ^b)\kb.i\\b^b) + 2iik]|Il2(b) 


— ' 2 






2 ' ' 


< C 


/2 ( ?,-l/2 
Kb,i 

2 




Vh 


L'^{Kb,i) 


Vpvvh 

h-^!^ 


Kb,\ 




V^v 


L^(Kb.i) 

Therefore, for any (5 > 0, we have 

(4.4) ^ [ {pVvh ■ ns} [vh]ds 

„ A, Jb 

B£SI 


LHKb, 2 ) 


L^{Kb,2 

1/2 


1 


1 


\B\ 


V2 


lHKb.2)J \b\ 

L^(B) ■ 


/2 llb/»]llL2(_B) 


- E ^ 

Best ^ 

<1 E I ' 


L'^{Kb,i) 


V^V 


Vh 


L'^{Kb, 2 ) 


1/2 


\B\ 


V2 


LHB) 


KeTh 


L^{K) 25 




Beet 


Then for e = 0 we let d = 1 and = C, and for e = —1 we let i5 = 1/2 and 
(7^ = 5(7/2, where C is in the above inequality. The coercivity result (4.2) follows 
from using these parameters in (4.4) and putting it in (4.3). □ 

In the error estimation for the IFF solution, we need to use the fact that both 
linear and bilinear IFE spaces have the optimal approximation capability [29, 30, 38]. 
In particular, for every u G satisfying the interface jump conditions (1.3) and 

(1.4), there exists a constant C such that the interpolation IhU in the (either linear 
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or bilinear) IFE space 5'^(n) has the following error bound: 

(4.5) \\u — Ihu\\j^ 2 (p^hi ||m —j < Ch ||w||^ 2 (f 2 ) ■ 

\TeTh / 

In addition, we also need the error bound for I^u on interface edges which is given in 
the following lemma. 

Lemma 4.2. For every u € satisfying the interface jump conditions (1.3) 

and (1.4), there exists a constant C independent of the interface such that its inter¬ 
polation IhU in the IFE space Sh{Il) has the following error bound: 

(4.6) ||/3(V(u - Ihu))\K ■ nB||i2(B) < C{h^ ||u|||3(o) + h \\ u \\% 2 (^ k ) ) 

where K is an interface element and B is one of its interface edge. 

Proof. We give a proof for linear IFEs, and the arguments can be used to establish 
this error bound for bilinear IFEs. 

Without loss of generality, let K = AA 1 A 2 A 3 be an interface triangle such that 

(4.7) 24i = (0,/r),242 = (0,0),243 = (/i,0) 
and assume that the interface points on the edge of K are 


(4.8) 


D = (0, d), E = {e,h — e) 


with Ai G K'^. Also we only discuss B = A 1 A 2 , the estimate on the other interface 
edge can be established similarly. 

By Lemma 3.3 and Lemma 3.4 in [38], for every X G DA 2 , we have 

dMx) 


{Ihu{X) - u{X))p = {N-{D) - Nj^)Xu-iX)iA, - D)- 


dp 


(4.9) 

where 


Xi[X )—^- \-l 2 {X )— -h/3(A)—j;-—, p = x,y, 


dp 


dp 


N-{D) = 


_ / '"'y 


{D)'^ + pna;{D)'^ {p-l)nx{D)ny{D) 

{p - l)nx{D)ny{D) ny^D)"^ + pn^iD)'^ 




2 ) 


dp 

+ pnl {p - l)h^hy 
{p - l)n^ny nl + pnl 


p = P /, n(X) = (nx{X),ny{X)y is the normal to F at X, n{DE) = {ux^Uy) 


is the normal of DE, and 


h{X) = il-td){N-iD)-I) 


dS/u 

dt 


{tD + (1 - t)X) ■ (Ai - X)dt 


d‘^ u~^ 

(4.10) +/ il-t)—^{tAi + {l-t)X)dt+ {l-t)—^itAi + {l-t)X)dt 

Jo dt Jt^ dt 


(4.11) h{X)= [ {l-t)^^{tA 2 + {l-t)X)dt 

Jo dP 


i = 2,3, 
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where D = tj^Ai + (1 — t^)X = X + ~ ^)- By Lemma 3.1 and Theorem 2.4 of 

[38], we have 

(4-12) j_{{N-{D)-NTy^)Vu-{X){A^-D)^^p-^ dX < |H|^3(n-), 

for p = x,y. By direct calculations we have 


dX/u 

dt 


{tD + (1 - t)X) ■ {Ai 


X) 


< 


u^^{tD + (1 — t)X){xd — x){xi — a;)| 


+ \ua:y{tD + (1 - t)X){yd - y){xi - a;)| + \uy^{tD + (1 - t)X){xd - x){yi - y)\ 
+ \u~y{tD + (1 - t)X){yd - y)(yi - y)|), 


and 

<( + (1 - i)X){xi - x){xi - a;)| + \uly{tA, + (1 - t)X)(yi - y){xi - x)\ 

+ + (1 - t)X)ixi - x){y^ - y)\ + \ulyitAi + (1 - t)X){yi - y){y, - y)|) 

where s = ±,i = 1,2,3. Let Iid{X),i = 1,2,3 be three integrals in Ii{X), respec¬ 
tively. Then, by Theorem 2.4 of [38], we have 

dX < ^ {1-tdf \uyy{OAyd + il-t)y){yd-y)ih-y)\^ dtdy 
(4.13) <Ch^ r\u-y{0,z)\^dz<Ch^\\u\\l,^^.^, 


dX<cj \uyy{G,y + t{h-y))\^{h-yf{l-tfdtdy 

(4.14) <Ch^ /'%-(0,z)|'dz<C'h2||u||^3(o-), 

J 0 


J__(^h,3iX)^^^^^^ dX<cJ \u^y{Q,y + t{h-y))\^ {h-yf{l-tfdtdy 


' da-2 
(4.15) 

Similarly, we can show that 


<Ch2 / \uyy{^x)\^ dz<Ch^\\u\\\^^+). 
J d 


dX<cJ^ J \uyy{0,{l-t)y)\'^ y‘^{l-tfdtdy 
< Ce (0, z)|" dz < Ch^ ||«||^ 3 (o-) ■ 


(4.16) 


^0 
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For the term involving / 3 (X), we have 


(4.17) 


hix) 


dMX)\ 


SAT V 9 p J 


dX 


<c(h'^[ [ \u^^{th,{l-t)y)\^ {l-t)'^dtdy+ [ [ \u^y{th,{l - t)y)\^ y'^{l - t)'^dtdy 

\ Jo Jo Jo Jo 

n l pd pi \ 

\uya;{th,{l-t)y)\^y'^{l-tfdtdy+ / / \u~y{th,{l - t)y)\^ y'^{l - t)'^dtdy] . 

Jo Jo j 


Let th = p, {1 — t)y = q, then we have 


and. 


Hence, 


P _ g _ g 

1 — t 1 — p/h’ 


ri — p fi — tri ri ri 


-dpdq 


J J \uxx(.th,{l-t)y)\^ {1-tfdtdy = jj \u^^{p,q)\^ ^ 

A_DA2 ^3 

< JJ \uxxiP,q)\^ {>^-p)dpdq<Ch\\u\\%2(^K)^ 

ADA2A3 

J^ J^ \uxy{th,{l-t)y)\^y^{l-tfdtdy= JJ \u-y{p,q)\^ j^^dpdq < Ch\\u\\%2^ii) 


/o ^0 

rd pi 


DA 2 A^ 


J^ J^ \'^vvitl^A^-t)y)\^y'^0-tfdtdy= JJ \uyy{p,q)\^ j^^dpdq<Ch\\u\\%2^K)- 


Using these estimates in (4.16), we have 

(4.18) dX < Ch\\u\\%2^j,y 

Finally, the inequality (4.6) follows from putting estimates (4.12)-(4.16), (4.17) into 
(4.9). □ 

Remark 4.1. For every u € satisfying the interface jump conditions 

(1.3) and (1.4), we can also show that there exists a constant C independent of the 
interface such that the interpolation IhU in the IFF space Sh{^) fulfills 

(4.19) ||/3(V(m — Ihu))\K ■ < Ch^ , 

where K is an interface element and B is one of its interface edge. A proof for (4.19) 
is given in Appendix A which uses arguments similar to those for proving Lemma 

4 . 2 . 
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Now, we are ready to derive the error bound for IFE solutions generated by the 
partially penalized IFE method (2.4). 

Theorem 4.3. Assume that the exaet solution u to the interface problem (1.1)- 
(1.4) is in H^{n) and Uh is its IFE solution generated with a = \ on a Cartesian 
(either triangular or rectangular) mesh 7h- Then there exists a constant C such that 

(4.20) \\u — Uh\\f^ < Ch\\u\\j^ 3 ^^y 


Proof. From the weak form (2.3) and the IFE equation (2.4) we have 

(4.21) ah{vh,Uh - Wh) = ah{vh,u - Wh), Vvh^Wh € 

Letting vt = ut — Wh in (4.21) and using the coercivity of aft(-, •), we have 

(4.22) K \\uh - WhWl ^ \ah{uh - Wh,Uh - Wh)\ = \ah{uh -Wh,u- Wh)\ 

/ l3V{uh-Wh)-V{u-Wh)dX + / {ldV{u-Wh)-'n.B}[uh-Wh]ds 

JK „ Jb 


< 


Ken 


Beet 


/ {(dV {uh - Wh) ■-O-b} [u - Wh]ds + y] / -^a[uh - Wh][u - Wh]ds 
We denote the four terms on the right of (4.22) by Qi, i = 1, 2, 3,4. Then, 

1/2 / s 1/2 


Qi < f y] I/3^^^V(u - w/t) ^ f y] \ld'^/'^V{uh-Wh) 

\Ken \Ken 

< ^max(^-,/3+)||V(u-w;?,)||^2(n) + f \ - Wh) 


Ken 


L^K) j 
2 

LHK) 


(4.23) < C ||V(u - w/.)||i2(o) -f - \\uh - WhWl ■ 


E ^IIK-^^]fL^(B)+C'E ^ll{/5V(u-u;,)-nB}||Es) 

Beii Beil ^ 

(4.24) < J llii, - WhWl + C E ^ ll{/3V(« - Wh) ■ nB}\\lyB) ■ 

BeSi ^ 

To bound Qsj for each B G we let KB,i G Th,i = 1,2 be such that B = Kb,iCKb, 2 - 
First, by the standard trace inequality on elements for functions, we have 

||[u-Wh]|li2(B) < \\{u - Wh)\KB,i\\n(B) + Wi'^ - '*^h)\KB, 2 \\n(B) 

< Ch ^||u — Wh\\B'i(^KB,i) fo ^ ~ 

+ (||m - WhWh^Knn + ^ “ '^^)llL=(ifB,2)) ■ 
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Then, applying the trace inequalities established in Lemma 3.2 or Lemma 3.5 depend¬ 
ing on whether linear IFEs or bilinear IFE are considered, we have 


\\{l3V{Uh-Wh)-Ta.B}\\L^(B) 

< ^ ^/PW{Uh- Wh) 


L^(Kb,i) 


^/pV{uh - Wh) 


l^(kb, 2 : 


Hence 


Qa < |e| \\{l3^{uh-Wh)-nB}\\B2^B)\\['^-'^h]\\L2(B) 


Beei 


(4.25) < Ch 2 ||V(u - ^ \\uh - Wh\\l ■ 

To bound Qi, by the standard trace inequality, we have 


f cr° 


LHB) 


< 


< 


'B 


\B\ 


\Br 


^||(m T II (m |iCB,2 II 

c \B\ \Kb,i\~^ ( ||u - Wh\\B2(^KB,i) + ^ ll^(“ “ '^h)\\L2(KB,i) Y 


+ T^C* \B\ \Kb,2\ M II« - ^h\\L^^KB,2) + ^ ll^(“ “ '^h)\\L^(^KB,2) Y 


< Ch {\\u - Wh\\B2(^KB,i) + - '^h)\\L2(KB,i)Y 

(4.26) + (||u - Whh2^KB.2) + h l|V(« - ^>^)\\L-i.KB,2) )'> 

where we have used the facts that 


\B\ = h or \B\=y/2h, \KB,i\ = -j^^ or \KB,i\ = = 1,2. 


Then 


Qa<Y. 

Beii 


K a 


B 


6/p B 


K, 1,0 3 

-6 "“‘-'"'■"‘+ 2 ; 


3 /* (T^ X 

[uh - Wh] [uh “ ^ J [■“ “ “ Wh]dsj 

r jjO 

X] J J^bj'- Wh][u - Wh]ds 


(4.27) < '^\\uh-Wh\\l + Ch ^°‘^^Y\\u-Wh\\l2^Q)+h^\\V{u-Wh)\\l2(^^^). 
Then, we put all these bounds for Qi, i = 1, 2, 3,4 in (4.22) to have 


\\uh - WhWl ^ C' ||V(u - Wh)\\l2^^^ + ‘^ XI ll{^^(« - '^h) ■ nB}||^2(B) 


BGSi 


+ Ch 2 (||u - Wh\\l2(^a) + ll^(« - w^'»)llL2(n)) 

+ Ch-<^°‘+^'> {\\u - Wh\\l 2 ^Q^ + h'^ l|V(u-'u;h)|li 2 (n)). 
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Hence, we let Wh = IhU in the above and using the optimal approximation capability 
of linear and bilinear IFE spaces (4.5) to have 


(4.28) 

\\uh - Ihu\\l < C{h? + ^ X! ll{/3V(« - Ihu) ■ nB}||i2(B) ■ 

For the second term on the right in (4.28), we use Lemma 4.2 to bound it: 

(4.29) < ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

Here we have used the fact that the number of interface elements is of 0{h~^). Hence, 
for a = 1, we can combine (4.28) and (4.29) to have 

(4.30) \\uh — Ihu\\f^ < Ch\\u\\jj3f^Q^ . 

In addition, using the optimal approximation capability of linear and bilinear IFE 
spaces (4.5) and (4.26), we can show that 

(4.31) \\u — Ihu\\f^ < Ch\\u\\fj2^p^y 

Einally, the error estimate (4.20) follows from applying (4.30) and (4.31) to the fol¬ 
lowing standard inequality: 


||m - UhWh < ll« - IhuWt^ -f \\uh - ■ 


□ 

Remark 4.2. If the exact solution u to the interface problem (1.1)-(1.4) is in 
the function space tF^’“(r2), then, using Remark 4-1 and arguments similar to those 
for the proof of Theorem 4-3, we can show that the IFE solution Uh generated with 
a = I on a Cartesian (either triangular or rectangular) mesh 7h has the following 
error estimate: 


(4.32) ||m R/tll/j < ||u||^2(q) + ). 


5. Numerical Examples. In this section, we present a couple of numerical 
examples to demonstrate features of the partially penalized lEE methods for elliptic 
interface problems. 

Eor comparison, the interface problem to be solved in the numerical examples is 
the same as the one in [30]. Specifically, we consider the interface problem (1.1)-(1.4) 
except that (1.2) is replaced by the non-homogeneous boundary condition u|an = <?• 
Let the solution domain H be the open rectangle (—1,1) x (—1,1) and the interface 
F be the circle centered at origin point with a radius tq = 7r/6.28 which separates fl 
into two sub-domains, denoted by CL~ and i.e., 

^~ = {{x,y) ■■ + y^ <rl}, and H+= {(a:, y) : -b > Tq}. 
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The exact solution u to the interface problem is chosen as follows 


(5.1) 


u{x,y) 




if r < To, 
otherwise, 


where a = 5, and r = + 2/^- The functions / and g in this interface problem 

are consequently determined by u. The Cartesian meshes Tinh > 0 are formed by 
partitioning fl into N x N congruent squares of size h = 2/N for a set of values of 
integer N. 

To describe our numerical results, we rewrite the bilinear form in the partially 
penalized IFE methods (2.4) as follows 


ahiuh,Vh) = ^ / pVuh-VvhdX + 5 ^ / {jSVuh ■'n.B}[vh]ds 

Ten 

(5.2) XJ y {ld'^Th-nB}[uh]ds+ XJ ^ bj-h][vh\ds. 


Best 


Best 


When (5 = 0, e = 0, and cr^ = 0 for all B G the partially penalized bilinear 
IFE method becomes the classic bilinear IFE method proposed in [29, 30]. When 
(5 = — 1, and > 0, we call the partially penalized IFE methods corresponding to e = 
— 1, 0,1, respectively, the symmetric, incomplete, and nonsymmetric PPIEE methods 
because of their similarity to the corresponding DG methods [50]. In our numerical 
experiment, we use the parameter a = 1 as suggested by our error estimation. Also, 
we use a^B = 10max(/3“,/3+) in the SPPIPE and IPPIPE schemes and = I in the 
NPPIFE scheme. 

In the first example, we test these lEE methods with the above interface problem 
whose diffusion coefficient has a large jump (/3“,,5+) = (1,10000). Because of the 
large value of the exact solution u{x,y) varies little in which might be one of 
the reasons this example is not overly difficult for all the IFE methods, and they indeed 
perform comparably well, see Tables 5.1, 5.2, and 5.3. Specifically, all the partially 
penalized IFE methods converge optimally in the i/^-norm as predicted by the error 
analysis in the previous section. The classic IFE method also converges optimally in 
the even though the related error bound has not been rigorously established 

yet. The data in Table 5.2 demonstrate that all the IFE methods converge in the 
L^-norm at the expected optimal rate. Even though they all seem to converge in the 
L°°-norm, the data in Table 5.3 do not reveal any definite rates for them. 


Table 5.1 

Comparison of \uti — for different IFE methods with /3~ = 1, /3+ = 10000. 



Classic IFE 

SPP IFE 

IPP IFE 

NPP IFE 

N 

l-lni 

rate 

\-\m 

rate 

\-\m 

rate 

l-lni 

rate 

20 

1.9495B-2 


1.9538E-2 


1.9538E-2 


1.9482E-2 


40 

1.0539E-2 

0.8873 

1.0601E-2 

0.8821 

1.0602E-2 

0.8820 

1.0562E-2 

0.8832 

80 

5.4219E-3 

0.9590 

5.5838E-3 

0.9249 

5.5844£;-3 

0.9248 

5.4945E-3 

0.9428 

160 

2.7157T;-3 

0.9975 

2.7720F;-3 

1.0103 

2.7724E-3 

1.0103 

2.7291E-3 

1.0096 

320 

1.3535E-3 

1.0046 

1.3887E-3 

0.9973 

1.3890E-3 

0.9970 

1.3613E-3 

1.0035 

640 

6.7876E-4 

0.9957 

6.9915E-4 

0.9900 

6.9923E-4 

0.9902 

6.7915E-4 

1.0032 

1280 

3.4023£;-4 

0.9964 

3.4439E;-4 

1.0216 

3,4447E-4 

1.0214 

3.3979E-4 

0.9991 

2560 

1.7079£;-4 

0.9942 

1.7036E;-4 

1.0155 

1.7038E-4 

1.0156 

1.6994E-4 

0.9996 
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Table 5.2 

Comparison of \\uf^ — Li||^2(fj) /o’" different IFE methods with fi~ = 1, /3+ = 10000. 



Classic IFE 

SPP IFE 

IPP IFE 

NPP IFE 

N 

II • IIl2 

rate 

II • IIl2 

rate 

II • IL 2 

rate 

II ■ IIl2 

rate 

20 

1.1175E-3 


1.1273E-3 


1.1276E-3 


1.1179E-3 


40 

2.8572£;-4 

1.9676 

2.9171E-4 

1.9503 

2.9181E-4 

1.9501 

2.8567E-4 

1.9683 

80 

7.5990B-5 

1.9107 

8.5150E-5 

1.7764 

8.5223E-5 

1.7757 

7.6342E-5 

1.9038 

160 

1.8116B-5 

2.0685 

2.2589E-5 

1.9144 

2.2645E-5 

1.9120 

1.8098E-5 

2.0766 

320 

4.4753E-6 

2.0172 

6.1332E-6 

1.8809 

6.1629E;-6 

1.8775 

4.5193E-6 

2.0017 

640 

1.0969T;-6 

2.0286 

1.6502F;-6 

1.8940 

1.6502E-6 

1.9010 

1.1235E-6 

2.0081 

1280 

2.6689E-7 

2.0391 

3.7104E-7 

2.1530 

3.7275E-7 

2.1464 

2.7680E-7 

2.0211 

2560 

6.3940E-8 

2.0615 

7.3251E-8 

2.3407 

7.3657E-8 

2.3393 

6.8809E-8 

2.0082 


Table 5.3 

Comparison of \\ufi — Li||Loo(f 2 ) for different IFE methods with f)~ = 1, /3+ = 10000. 



Classic IFE 

SPP IFE 

IPP IFE 

NPP IFE 

N 

II ■ Hi/®® 

rate 

II ■ III’®® 

rate 

II ■ lU- 

rate 

II-lU- 

rate 

20 

8.8830E-4 


9.9803E-4 


9.9820E-4 


8.4663E-4 


40 

4.3525E-4 

1.0292 

5.2256E-4 

0.9335 

5.2185E-4 

0.9357 

4.2059E-4 

1.0093 

80 

1.6536F;-4 

1.3962 

2.2692F;-4 

1.2034 

2.2745E-4 

1.1981 

1.6816E-4 

1.3225 

160 

7.4603F;-5 

1.1483 

9.8441F;-5 

1.2049 

9.8619E-5 

1.2056 

7.2067E-5 

1.2225 

320 

1.2972F;-5 

2.5238 

3.6045F;-5 

1.4494 

3.6421E-5 

1.4371 

1.4362E-5 

2.3271 

640 

6.1332E-6 

1.0807 

1.7607E-5 

1.0337 

1.7531E-5 

1.0548 

6.4702E-5 

1.1504 

1280 

2.5136E;-6 

1.2869 

5.1961E;-6 

1.7606 

5.1799E-6 

1.7589 

1.6993E-6 

1.9289 

2560 

1.2810E;-6 

0.9725 

1.0242£;-6 

2.3430 

1.0274E-6 

2.3339 

9.1534E-7 

0.8926 


However, we have observed that the partially penalized IFE methods outperform 
the classic IFE method in many situations. We demonstrate this by numerical results 
generated by these IFE methods for the interface problem whose diffusion coefficient 
has a typical moderate jump (/3“,/3+) = (1,10). IFE solution errors are listed in 
Tables 5.4, 5.5, and 5.6. From the data in Table 5.4, we can see that all the partially 
penalized IFE methods maintain their predicted 0{h) convergence rate in the H^- 
norm over all the meshes up to the finest one, while the classic IFE method slightly 
looses its convergence rate in the i/^-norm when the mesh becomes very fine. The 
effects of the penalization are more prominent when the errors are gauged in the L^- 
norm and T°°-norm. According to the data in Table 5.5, all the partially penalized 
IFE methods converge at the optimal rate in the L^-norm but the L^-norm rate of 
the classic IFE method clearly degenerates when the mesh becomes finer. Similar 
phenomenon for the L°“-norm convergence can be observed from those data in Table 
5.6. 


Table 5.4 

Comparison of \uh — for different IFE methods with fi~ = 1, /3+ = 10. 



Classic IFE 

SPP IFE 

IPP IFE 

NPP IFE 

N 

l-lni 

rate 

l-lm 

rate 

l-lm 

rate 

\-\hi 

rate 

20 

6.4758E-2 


6.4751E;-2 


6.4753E-2 


6.4719E-2 


40 

3.2779E-2 

0.9823 

3.2650F;-2 

0.9878 

3.2651E-2 

0.9878 

3.2637E-2 

0.9877 

80 

1.6596E-2 

0.9723 

1.6386E-2 

0.9947 

1.6386E-2 

0.9947 

1.6382E-2 

0.9944 

160 

8.4072F;-3 

0.9811 

8.2081F;-3 

0.9974 

8.2082E-3 

0.9973 

8.2063E-3 

0.9973 

320 

4.2566F;-3 

0.9819 

4.1071F;-3 

0.9988 

4.1071E-3 

0.9989 

4.1068E-3 

0.9987 

640 

2.2267E-3 

0.9348 

2.0544E-3 

0.9994 

2.0545E-3 

0.9994 

2.0544E-3 

0.9994 

1280 

1.1795E-3 

0.9167 

1.0274£;-3 

0.9997 

1.0274E-3 

0.9997 

1.0274E-3 

0.9996 

2560 

6.6893E-4 

0.8182 

5.1377E-4 

0.9998 

5.1377E-4 

0.9998 

5.1377E-4 

0.9998 


It is known that the point-wise accuracy of the classic IFE methods in the liter- 
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Table 5.5 

Comparison of ||ji^ — for different IFE methods with fi~ = 1, /3+ = 10. 



Classic IFE 

SPP IFE 

IPP IFE 

NPP IFE 

N 

II • IIl2 

rate 

II • 111,2 

rate 

II • IL 2 

rate 

II ■ IIl2 

rate 

20 

4.3003E-3 


4.2945E-3 


4.2989E-3 


4.2869E-3 


40 

1.0622£;-3 

2.0174 

1.0749E-3 

1.9983 

1.0745E-3 

2.0003 

1.0626E-3 

2.0124 

80 

2.6196B-4 

2.0196 

2.6833E-4 

2.0021 

2.6797E-4 

2.0036 

2.6440E-4 

2.0067 

160 

6.4952£;-5 

2.0119 

6.7047E-5 

2.0008 

6.6872E-5 

2.0026 

6.5876E-5 

2.0049 

320 

1.6311E-5 

1.9935 

1.6829E-5 

1.9942 

1.6794£;-5 

1.9935 

1.6594E-5 

1.9891 

640 

4.4482B-6 

1.8746 

4.2038F;-6 

2.0012 

4.1934E-6 

2.0017 

4.1383E-6 

2.0035 

1280 

1.4445E-6 

1.6226 

1.0501E-6 

2.0012 

1.0472E-6 

2.0015 

1.0336E-6 

2.0014 

2560 

6.7593E-7 

1.0956 

2.6254E-7 

1.9999 

2.6149E-7 

2.0017 

2.5821E-7 

2.0357 


Table 5.6 

Comparison of \\ufi — for different IFE methods with l3~ = 1, = 10. 



Classic IFE 

SPP IFE 

IPP IFE 

NPP IFE 

N 

II ■ \\l°° 

rate 

II ■ Hi-®® 

rate 

II ■ IIl°° 

rate 

II • IIl°° 

rate 

20 

1.0969E-3 


1.3680E-3 


1.3785E-3 


1.0082E-3 


40 

5.4748F;-4 

1.0026 

3.9775F;-4 

1.7822 

3.9769E-4 

1.7934 

1.9172E-4 

2.3947 

80 

5.0812E-4 

0.1077 

1.0601E-4 

1.9077 

1.0582E-4 

1.9100 

5.4491E-5 

1.8149 

160 

2.2635E;-4 

1.1667 

3.1598E-5 

1.7463 

3.1217E-5 

1.7612 

1.4045E-5 

1.9559 

320 

1.2290E;-4 

0.8811 

7.0324£;-6 

2.1677 

6.9364E-6 

2.1701 

3.5092E-6 

2.0009 

640 

7.0810E-5 

0.7954 

1.9288E-6 

1.8664 

1.9213E;-6 

1.8521 

9.1942E-7 

1.9324 

1280 

3.4111E-5 

1.0537 

5.0505E-7 

1.9332 

4.9869E;-7 

1.9459 

2.2932E-7 

2.0034 

2560 

1.7815E-5 

0.9371 

1.0199E-7 

2.3080 

1.0075E-7 

2.3074 

5.9784E-8 

1.9395 


ature is usually quite poor around the interface and we suspect the discontinuity of 
the IFE functions across interface edges is the main cause for this shortcoming. With 
the penalty to control the discontinuity in IFE functions, a partially penalized IFE 
method has the potential to produce better point-wise approximations. To demon¬ 
strate this, we plot errors of a classic bilinear IFE solution and a NPP IFE solution in 
Figure 5.1. The IFE solutions in these plots are generated on the mesh with 80 x 80 
rectangles. From the plot on the left, we can easily see that the classic IFE solution 
has much larger errors in the vicinity of the interface. The plot on the right shows 
that the magnitude of the error in the NPP IFE solution is much smaller uniformly 
over the whole solution domain. This advantage is also observed for other partially 
penalized IFE solutions. 



Fig. 5.1. A comparison of point-wise errors of a classic bilinear IFE solution and a nonsym- 
metric partially penalized IFE solution. The coefficient in the interface problem has a moderate 
jump: {0-,^+) = (1,10). 
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A. A proof of Remark 4.1. We give a proof for the linear IFEs, and the same 
arguments can be used to show this error bound for bilinear IFEs. 

Without loss of generality, let K = AA 1 A 2 A 3 be an interface triangle whose 
vertices and interface intersection points are given in (4.7) and (4.8) with Ai G iF+. 
Also we only discuss B = A 1 A 2 , the estimate on the other interface edge can be 
established similarly. 

By Lemma 3.3 and Lemma 3.4 in [38], for every X G ZIA 2 , we have 


{Ihu{X) 

(A.l) 


i(A))p = {N-{D) - iV^)Vu-(A)(Ai - D) 


d<^i{X) 

dp 


■ hix ) 


dMX) 

dp 


-hiX) 


dMX) 

dp 


■MX) 


dMX) 

dp ’ 


p = x,y, 


where N~{D), N-^, p, n(A), and n{DE) are defined the same as in the proof of 
Lemma 4.2. The quantities Ii{X), i = 1,2,3, are given in (4.10) and (4.11). By 
Lemma 3.1 and Theorem 2.4 of [38], we have 


(A.2) 


J_ (^{N-{D) - iV^)Vu-(A)(Ai - dX < Ch^ Ijulj"^ 




for p = x,y. By direct calculations we have for i = 1,2,3,p = x, y 
dX/u~ 


dt 

d'^u^ 


-{tD+{l-t)X)-{Ai-X) 


< ( I'^xxi'tD + (1 — t)X) \ + \uMidd + (1 — t)X) 


dE 


[tAi + (1 — t)X) 


+ + (1 — t)X)\ + \uyy{tD + (1 — t)A)| )h?, 

^ + (1 ~ ^)X)\ + \u^y{tAi + (1 — t)X)\ 

“b (1 ~ “b + (1 ~ ^)^)| ) ■ 


By these estimates and Theorem 2.4 of [38], we have 


(A.3) 


MX) 


dMX) 


AWT V dp 


dX < Ch? ||wj|^y2,oo(Q) , 


Then, using (A.2) and applying (A.2) and (A.3) we have 

(A.4) [_ _((4m(A) - u(A))p)^(iA < 11^11^2, oo(n) , p = x,y. 

J DA2 

Similar arguments can be used to show 

(A.5) f _ {{Ihu{X)-u{X))p)'^dX <CM\\u\\'^ 2 ,^,q), P = x,y. 

JAiD 


Einally, the estimate (4.19) on the interface edge B = A 1 A 2 follows from (A.4) and 
(A.5). 
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